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UNITED STATES SPECIFTCATTOKT 
TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN THAT WE, EMANUELA KELLER, ANDREAS NADLER and 
PETER NIEDERER, all citizens of Switzerland, residing at 
Hornhaldenstrasse 40, 8802 Kilchberg, Switzerland, 
Eintrachtstrasse 11, 8820 Wadenswil, Switzerland, and 
Nordstrasse 366, 8037 Zurich, Switzerland, respectively, have 
invented certain new and useful improvements in a 

DEVICE FOR MEASURING THE BLOOD FLOW IN AN ORGAN 

of which the following is a specification. 
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BACKGROUND OF THE INVENTION 



CROSS REFERENCE TO RELATED APPLICATIONS 

Applicants claim priority under 35 U.S.C. §119 of 
German Application No. 103 15 574.0 filed April 5, 2003. 

1. Field of the Invention 

The present invention relates to a device for measuring 
the blood flow in an organ, using an injected indicator. 

2 . The Prior Art 

Devices for measuring the blood flow in an organ, 
particularly the cerebral blood flow (CBF) , using an 
injected, essentially inert indicator, have been known for a 
long time, and some of them are in clinical use. Devices 
based on older methods, such as the xenon dilution technique 
(Obrist, W.D., Thompson, H. K. , Wang, H.S., and Wilkinson, 
W.F. 1975, Regional cerebral blood flow estimated by 133 
xenon inhalation. Stroke 6: 245-256), are often difficult to 
implement, in technical terms, and prove to be time-consuming 
to use. Recently, technologies that use near-infrared 



spectroscopy (NIRS) in combination with indocyaningreen (ICG) 
as an indicator, in order to investigate or monitor cerebral 
blood flow, have gained in importance. 

Generally, the use of such devices requires the invasive 
measurement of the input function, i.e. the arterial 
concentration of the indicator over time. In this 
connection, arterial fiber optic catheters are used, for 
example. 

Most recently, the use of non-invasive measurement 
methods for monitoring cerebral blood flow have also been 
proposed (Keller, E. , Wolf, M. , Martin, M. , Fandino, J. and 
Yonekawa, Y. 2001. Estimation of cerebral oxygenation and 
hemodynamics in cerebral vasospasm using indocyaningreen 
(ICG) dye dilution and near infrared spectroscopy (NIRS) . A 
case report. J. Neurosurg . Anesthesiol . 13: 43-48). In this 
connection, optodes are affixed to the head, which function 
as transmitters and receivers of the near infrared radiation, 
by means of which the indicator concentration in the cerebral 
vascular system is determined. The algorithms used for 
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evaluation are described in the literature (Keller, E. , 
Nadler, A., Imhot, H.-G., Niederer, P., Roth, P. and 
Yonekawa, Y. 2002. New Methods of Monitoring Cerebral 
Oxygenation and Hemodynamics in Patients with Subarachnoid 
Hemorrhage. Acta Neurochir. [Suppl.] 82: 87-92). 

A significant problem of conventional NIRS devices is 
that their evaluation units work with a certain inaccuracy, 
which arises because in the evaluation, until now, the point 
of departure has been the simplifying assumption that no flow 
of the indicator out of the system being investigated takes 
place during the first rise in the input signals (the 
so-called "rise time") . The accuracy of the results 
delivered within the scope of the evaluation therefore 
decisively depends on the sharpness of the first peak of the 
measurement signal after administration of the indicator. 
Since injection of the indicator cannot be performed in a 
time that is as short as desired, and at a location that is 
as close to the measurement location as desired, for obvious 
practical reasons, the dependence of the quality of the 
measurement results on the sharpness of the first peak of the 
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measurement signal after administration of the indicator was 
accepted as being fundamental, until now. 

Alternative measurement methods, such as positron 
emission tomography (PET), single-photo photon emission 
computer tomography (SPECT) , or per fusion-weighted magnetic 
resonance spectroscopy can be implemented technically only by 
means of extremely expensive devices, and require transport 
of the patient in question to the measurement unit, in order 
to be used, and this transport often poses significant risk 
to the patient; in other words, the devices mentioned cannot 
be used as so-called "bedside" devices at the patient's bed. 

SUMMARY OF THE INVENTION 
It is an object of the present invention to provide a 
device that can achieve greater accuracy in the determination 
of the blood flow in an organ, particularly the cerebral 
blood flow. In particular, the accuracy of the evaluation 
results should be less dependent on the general conditions, 
i.e. particularly on how the indicator injection was 
performed. Furthermore, the device should be non-invasive, 
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and also implementable with a reasonable technical 
effort/expenditure . 

According to one aspect of the present invention, this 
object is achieved with a device for measuring the blood flow 
in an organ, using an injected indicator which includes a 
radiation source for emitting near infrared radiation into 
tissue of te organ at a first location, a sensor for 
detecting a proportion of the emitted near infrared radiation 
that exits from the organ at a second location, and an 
evaluation unit that detects the proportion of emitted near 
infrared radiation that exits from the tissue of the organ as 
an input signal. The input signal is composed of a pulsatile 
component and a non-pulsatile component. The evaluation unit 
is set up, in terms of program technology, to perform the 
following steps: 

(a) determination of the concentration of injected 
indicator, with reference to the tissue of the organ, from 
the non-pulsatile component of the input signal, 
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(b) iterative determination of an inflow function i(t) 
that characterizes the blood flow through the organ, by 
variation of a mean transit time mtt until a stop criterion 
is reached, 

(c) determination of the concentration of injected 
indicator, with reference to the blood volume in the organ, 
from the pulsatile component of the input signal and the 
iteratively determined inflow function i(t), 

(d) calculation of the blood volume in the organ as a 
quotient of the concentration of injected indicator with 
reference to the tissue of the organ and the concentration of 
injected indicator with reference to the blood volume in the 
organ , and 

(e) calculation of the blood flow in the organ as a 
quotient of the blood volume in the organ and the mean 
transit time mtt when the stop criterion has been reached. 
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Advantageous embodiments of the device according to the 
invention are described below. 

One of the advantages of the present invention is 
that the sharpness of the first peak of the measurement 
signal after administration of the indicator is no longer 
decisive for the accuracy of the results achieved, something 
that is completely surprising to a person skilled in the art, 
since the evaluation unit also takes the flow of indicator 
out of the organ during the first rise in the input signal 
into consideration, because of its program technology set-up 
according to the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Other objects and features of the present invention will 
become apparent from the following detailed description 
considered in connection with the accompanying drawings. It 
should be understood, however, that the drawings are designed 
for the purpose of illustration only and not as a definition 
of the limits of the invention. 
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In the drawings, 

FIG. 1 shows schematically, the most important 
components of an embodiment for determining the cerebral 
blood flow of a patient, in accordance with the invention, 
and 

FIG. 2 shows a flow chart of the sequence for the 
determination of the cerebral blood flow from the input 
signal, implemented in the evaluation unit in terms of 
program technology. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
The device according to the invention shown 
schematically in FIG. 1 serves to determine the cerebral 
blood flow of a patient, for example an intensive-care 
patient in neurosurgery. A first optode 2 and a second 
optode 3 are attached to the head 1 of the patient, by means 
of an elastic band (not shown) , at an optimized distance from 
one another. A radiation source (not shown) for emitting 
near infrared radiation into the cerebral tissue of the 
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patient is arranged either in first optode 2 itself, or 
separately from it, for example in a common housing with the 
evaluation unit 4. Where the radiation source is separate 
from first optode 2, the near infrared radiation is passed by 
means of a light guide to first optode 2, where the radiation 
is emitted. The wavelength of the emitted radiation and the 
indicator used must be coordinated with one another. For the 
indicator usually used, indocyaningreen, a wavelength of 
about 805 nm (but definitely in the range between 7 80 and 
910 nm) is ideal. The intensity of the proportion of the 
infrared radiation that exits from the cerebral tissue at the 
location of second optode 3 is detected by second optode 3 
and passed to the evaluation unit 4. 

This proportion of radiation is higher or lower, 
depending on the concentration of the indicator injected into 
the cerebral tissue of the patient. The indicator 
concentration with reference to the cerebral tissue depends 
on the indicator concentration in the blood that flows 
through the tissue, as well as on the amount of the blood 
that flows through the tissue. The amount of blood that 
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flows through the tissue changes periodically, with the 
heartbeat, and for this reason the intensity signal received 
has a pulsatile component, i.e. a periodically varying 
component. The pulsatile component of the intensity signal 
has a non-pulsatile component superimposed on it. 

The indicator concentration in the blood that flows 
through the tissue changes over time, because the blood that 
flows out of the cerebral tissue has a different 
concentration from the blood that is flowing in. The 
distribution kinetics of the indicator in the cerebral 
vascular system are partly responsible for this difference in 
concentration. The indicator concentration in the 
recirculated blood (blood newly flowing in) is determined by 
the distribution kinetics of the indicator in the entire 
circulatory system and by the decomposition of the indicator 
in the liver. 

Evaluation unit 4 is set up, in terms of program 
technology, for performing the evaluation steps listed in 
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FIG. 2 in the form of a flow chart, which will be explained 
below. 

The optical density OD is formed from the intensity 
signal as a negative decadic logarithm of the transmission. 
In this connection, transmission is understood to be the 
quotient of the intensity of the detected near infrared 
radiation and the intensity of the emitted near infrared 
radiation (and not in the strictly physical sense, since the 
detected radiation includes scattered and reflected 
components) . 

The (time -dependent) optical density, which is 
essentially proportional to the indicator concentration in 
the tissue, is divided up into its pulsatile component and 
its non-pulsatile component. 

An inflow function i(t) that corresponds to the arterial 
input function of the brain is iteratively determined from 
the non-pulsatile component. The sequence of the iteration 
is shown in the left branch of the schematic shown in FIG. 2. 
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A suitable mean transit time mtt, for example mtt = 7 s, is 
chosen as a starting parameter. The mean transit time mtt is 
sometimes also referred to as the "pass -through" time, and is 
a characteristic dwell time that corresponds to the time that 
a volume element needs, on the average, in order to pass 
through the system being considered. 

In an initialization step, the time variable t, as well 
as the inflow function i(t) for the range t < 0, and the 
related outflow function o(t), are set to zero., 

The inflow function i(t) describes the proportion of the 
change in the concentration of the indicator in the cerebral 
tissue that comes from the amount of inflowing blood; the 
outflow function o(t) describes the proportion of the change 
in the concentration of the indicator in the cerebral tissue 
that comes from the amount of outflowing blood. The time 
immediately after the injection of indicator is t = 0. 

Each iteration step includes a step-by-step calculation 
of an approximation of the inflow function i(t) and an 
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approximation of the outflow function o(t), as well as the 
calculation of an approximation of the transport function 
g(t), so that the m th approximation of the inflow function 
i(t), the outflow function o(t), and the transport function 
g(t) is calculated with the mth iteration step (let m be the 
counting variable) . In a computing step, the value of the 
inflow function for the time t is calculated according to the 
balance equation 

i(t) = d/dt (C tissue (t)) + o(t - t k ) . 

In this equation, t k is a constant small time interval, so 
that the value of the outflow function at the time t - t k is 
to be inserted for o(t - t k ) . The term d/dt (C tissue (t)) 
expresses the change in the indicator concentration with 
reference to the cerebral tissue. 

Within the framework of the iterative calculation of the 
inflow function, only relative function values are required, 
at first, because of scaling that takes place later by means 
of absolute values of the pulsatile signal component. For 
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this reason, the optical density OD(t) is used for the 
concentration C tissue (t). Furthermore, the term d/dt 
(C tissue (t)) is linearized, so that the above calculation step 
is implemented in the form 

i(t) = OD(t) - OD(t - t k ) + o(t - t k ) . 

In the next step, the convolution integral 
o(t)=i (t) *g(t) with the transport function g(t) is used for 
calculating the value of the outflow function at the time t. 
The transport function is formed according to one of the 
usual equations, for example the one published in Hoeft, A., 
Schom, B. , Weyland, A., Scholz, M. , Buhre, W. , Stepanek, E. , 
Allen, S.J., and Sonntag, H. 1991. Bedside assessment of 
intravascular volume status in patients undergoing coronary 
bypass surgery/- Anesthesiology 81: 76-86, and depends on the 
mean transit time mtt. A suitable equation is* ■ 

t rmtt £ I 

g(t) = -=J .. e 

V2jia.t 



- 14 - 



In this equation, a is a constant parameter selected for the 
system on the basis of empirical values (which fundamentally 
describes the width of an assumed dwell time distribution) . 

In a next step, the time variable t is increased by the 
increment t k . If t is smaller than an end value t 2 , the loop 
is run through again, starting with the calculation step 

i(t) = OD(t) - OD(t - t k ) + o(t - t k ) . 

If, on the other hand, the end value has exceeded t 2 , a check 
is performed in the next step to see whether the function 
progressions of the inflow function i(t) and the outflow 
function o(t) are plausible. A plausibility criterion (i.e. 
a stop criterion of the iteration) may be that neither the 
inflow function i(t) nor the outflow function o(t) have 
values below a threshold value. It is suitable to select 
this threshold value to be greater than or equal to 0. 
Furthermore, the plausibility criterion can include the 
requirement that the inflow function i(t) can be represented 
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as the sum of a finite number of functions that are similar 
to the form of the transport function g(t) . 

If the function progressions of the inflow function i(t) 
and the outflow function o(t) are not plausible, the mean 
transit time mtt is adapted by an increment by means of a 
suitable increase or decrease, and another iteration step is 
performed, which again begins with the initialization step 
and proceeds as explained. 

If, on the other hand, the function progressions of the 
inflow function i(t) and the outflow function o(t) are 
plausible, the iteration is stopped. The method then 
continues with the function progression of i(t) that was 
received and the value of the mean transit time mtt that was 
received, as described below. 

The indicator concentration in the tissue is also 
determined from the non-pulsatile component of the time 
progression of the optical density OD, as shown in the right 
branch of the schematic shown in Fig. 2. For this purpose, 
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the non-pulsatile component of the time progression of the 
optical density OD is simulated in an interval t ± > 0 to t 2/ 
by regression by means of an exponential function, and the 
density is extrapolated back to the time t = 0 . The time t 2 
is selected so that complete mixing of the indicator with the 
blood is guaranteed, i.e. so that concentration peaks due to 
recirculation can no longer be determined. Furthermore, the 
value of the non-pulsatile component is determined in the 
range t < 0, i.e. OD(t<0). 

In a next step, the actual calculation of the indicator 
concentration immediately after injection of the indicator, 
with reference to the tissue, C tissue , takes place, according 
to the formula 

^tissue — 

[OD(t->0) - OD(t<0)] / alCG 

whereby alCG is the absorption coefficient of the indicator, 
and OD(t->0) corresponds to the value of the optical density 
OD from the function progression for t approaching 0 that has 
been extrapolated back. 
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Evaluation steps implemented in order to take into 
consideration the pulsatile component are listed in the 
center branch of FIG. 2. An envelope is adapted to the 
pulsatile component of the time progression of the optical 
density OD. Furthermore, the amplitude A OD (t<0) / in 
other words the amplitude of the pulsatile component of the 
optical density, is determined for t<0. 

With the absolute values obtained from the pulsatile 
signal component, from a range in which the amplitude is 
great and therefore less susceptible to noise (for example 
for t from 0 to 60 s) , the inflow function i(t) that has 
previously been determined by iteration is scaled, for 
example by minimizing the sum of the difference squares 
[i(t) - Aq D ( t ) ] 2 . In the final analysis, this scaling results 
in greater accuracy, because while values obtained from 
direct measurement for greater t values are increasingly 
distorted by noise, due to lower indicator concentrations, 
i(t) describes the decrease in the indicator concentration in 
the blood well, even for great t values. 
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Proceeding from a "time interval of t la > 0 to t 2a , in 
which the indicator has already been distributed well in the 
system (i.e. at sufficiently great t la and t 2a ) , the scaled 
inflow function i(t) is extrapolated back to the time point 
t = 0 (for example by means of an exponential function 
simulated by means of regression, in the time interval 
t la > 0 to t 2 ) . 

The indicator concentration in the blood, C blood , required 
for calculating the cerebral blood volume CBV is calculated 
according to the following formula, wherein Aq D ( t > 0) 
corresponds to the value of the amplitude of the optical 
density for t approaching 0 that was obtained from the 
function extrapolated back (which is less dependent on the 
sharpness of the first signal peak after administration of 
the indicator, as compared with a value determined by direct 
measurement) : 

c [ AppO^O) } ) <i Hb C Hb 
Bl00d lA OD (t<0) J a ICG • 
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In this equation, alCG is the absorption coefficient of 
the indicator, is the absorption coefficient of the 
hemoglobin, and is the hemoglobin concentration in the 
blood. 

The cerebral blood volume CBV is calculated as the 
quotient of an indicator concentration in the blood, C blood , 
and the indicator concentration with reference to the 
cerebral tissue, C tissue , i.e. 

CBV = C tissue / C blood . 

The cerebral blood flow CBF is determined as a quotient 
of the cerebral blood volume CBV and the mean transit time 
mtt # in other words according to the formula 

CBF = CBV / mtt. 

Accordingly, while at least one embodiment of the 
present invention have been shown and described, it is to be 
understood that many changes and modifications may be made 
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thereunto without departing from the spirit and scope of the 
invention as defined in the appended claims. 
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